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Many theories of physics beyond the standard model 
(SM) predict additional vector bosons (e.g. Z'). These 
include (but are not limited to) extended gauge theories 
(e.g. SO(IO)), Kaluza-Klein states of the gluon or of Z 
bosons, axigluons, and topcolor PHTU]. Previous analyses 
of Tevatron data have excluded such narrow resonances 
with masses less than 725 GeV/c^ [TTUn] . 

In this Letter we describe a search for narrow resonant 
states decaying to top-antitop pairs in proton-antiproton 
collisions using the CDF II detector at the Tevatron at 
Fermilab. We discuss the experimental signature of reso- 
nant tt production and how we are to distinguish it from 
standard model tt production. We describe the recon- 
struction technique used in this analysis and the statis- 
tical tests we perform in examining the data for any sign 
of the hypothetical resonant production, and we describe 
the systematic uncertainties related to this analysis. Ob- 
serving no evidence for resonant production of tt pairs, 
we set upper limits on the cross section times branching 
ratio for resonant tt production. 

At the Tevatron, the tt pair-production cross section 
has been measured with great precision: cr -^^j = 7.70 ± 
0.52 pb [12]. However, this degree of precision leaves 
open the possibility that non-SM physics gives rise to a 



fraction of the total tt production. 

We search for a heavy vector boson decaying to ti in 
the final state where one top quark decays semileptoni- 
cally {t—i'£i'b) and the other hadronically (t—i'qq'b) [TB] 
by examining the tt invariant mass spectrum of candi- 
date events, where the event kinematics have been re- 
constructed by applying the SM QCD "Matrix Element" 
for ti production and decay ^13j . The observed spectrum 
is then compared to templates - models of signal (i.e., 
Z' ^ti) and background processes (e.g. SM ti, W+i\ets, 
WW/WZ) in an unbinned maximum-likelihood fit. By 
fitting the data to these models, we extract upper limits 
on the Z' — )■ ti cross section times branching ratio. We 
consider scenarios where the Z' width is 1.2% of the pole 
mass (this has been the benchmark scenario for narrow- 
resonance searches in ti enriched samples [5|). 

The CDF detector is a general purpose, azimuthally 
and forward-backward symmetric multipurpose collider 
detector. A detailed description can be found in Ref. [T7]; 
here we summarize details of detector components im- 
portant to this analysis. The transverse momenta (pt) 
and track parameters of charged particles are measured 
by an eight-layer silicon strip detector [TSl - BD] and a 96- 
layer drift chamber (COT) [3T], both within a 1.4 T mag- 
netic field. The COT provides tracking coverage with 
high efficiency for \r]\ < 1 [22]. Electromagnetic [23] and 
hadronic |24| calorimeters surround the tracking system. 
They are segmented in a projective tower geometry and 
measure the energies of charged and neutral particles in 
the central region {\ri\ < 1.1). A plug tile calorimeter cov- 
ers the forward region (1.1 < \ri\ < 3.6). Each calorime- 
ter has an electromagnetic shower profile detector posi- 
tioned at the shower maximum. The calorimeters are 
surrounded by muon drift chambers [25]. The central 
muon detectors comprise four layers of drift chambers 
which cover the region \r]\ < 0.6. Forward muons, with 
0.6 < \i]\ < 1.0, are detected by an additional four layers 
of drift chambers. Gas Cerenkov counters |26j measure 
the average number of inelastic pp collisions per beam- 
crossing and thereby determine the beam luminosity. 

We select ti candidates in data corresponding to an in- 
tegrated luminosity of 4.8 fb^^ in the Icpton-fjets chan- 
nel [27] by requiring one isolated charged lepton — an 
electron (reconstructed using the central electromagnetic 
calorimeter) or a muon (reconstructed with the central 
or forward muon detectors). 

Primary leptons must have rapidity |r/| < 1. Electrons 
must have transverse energy Et > 20 GeV and muons 
must have transverse momentum px > 20 GeV/c. Four 
or more central jets {\ri\ < 2) with Et > 20 GeV are also 
required. One of these jets must have a secondary vertex 
displaced from the primary event vertex such that it is 
"tagged" by the SECVTX algorithm |28j as being consis- 
tent with the decay of a long-lived 6-hadron. Jets with 
6-tags are restricted to \r]\ < 1. Large missing trans- 
verse energy [53] is also required, ^t > 20 GeV. Jets 
corrected |30j for multiple pp interactions in the event, 
non-uniformities in the calorimeter response along rj and 



component 


4 jets 


> 5 jets 


non-W^ 


46.1 ±35.7 


15.7 ±12.2 


Z+light flavor 


6.4 ±0.5 


1.6 ±0.1 


W+light flavor 


32.9 ±8.5 


7.4 ±3.1 


Wbb 


51.5 ±12.6 


12.4 ±3.7 


Wcc 


27.7 ±6.6 


7.3 ±2.1 


Wcj 


14.0 ±3.3 


3.0 ±0.9 


single top 


8.9 ±0.4 


1.4 ±0.0 


diboson 


9.1 ±0.6 


2.4 ±0.1 


total non-ff 


196.6 ±39.5 


51.2 ±13.3 


SM it 


667.1 ±61.8 


225.2 ±21.0 



TABLE L Estimate of sample composition in terms of signal 
and background. 





central e* 


central /i* 


forward fi^ 


4 jets, 1 tag 


480 


221 


131 


4 jets, > 2 tags 


110 


33 


21 


> 5 jets, 1 tag 


164 


81 


41 


> 5 jets, > 2 tags 


47 


21 


16 



TABLE II. Number of events observed in each subsample for 
data corresponding to an integrated luminosity of 4.8 fb^^. 



any non-linearity and energy loss in the uninstrumented 
regions of the calorimeters. The ^t is corrected for the 
primary vertex location and for any high-px muons in 
the event. 

The contribution of non-tt backgrounds satisfying this 
selection has been derived in precision measurements of 
the ti production cross section [37], these are tabulated 
in Table |l] Also shown is the expected SM ti content, 
given the event selection. However, in this analysis the 
SM ti content is estimated as the difference between the 
number of observed events and the estimate of resonant ti 
events plus non-tt background events. Tabic \n\ indicates 
the number of events observed in each of the subsamples 
considered in this analysis. 

For each candidate event, we apply the ti hypothesis 
- the observed event kinematics are mapped to the par- 
ton level using the information available in the SM QCD 
"Matrix Element" for ti production and decay [13]. By 
constraining the event kinematics to the hypothetical ma- 
trix element, we can reconstruct with enhanced precision 
any kinematic quantity which can be expressed in terms 
of the parton momenta. In this analysis, we reconstruct 
the ti invariant mass probability density function (pdf) 
for each event. 

The probability for ti production and decay, given p, 
the four-momenta of the decay products of the ti system, 
is 



7r(p|mt) 



cr{mt) 



dZadZbfk{Za)fl{zb)dakl{p\mt,Za,Zb), 

(1) 
where fk{za) and fi{zb) are the parton density func- 



tions for incoming gluons or quarks with flavor k and 
I, with z-components of parton momenta, Za and Zb- 
daki{p\mt,Za,Zb) is the diff'erential cross section - it is 
proportional to the squared SM matrix element for ti 
production and decay. In this analysis, we fix the top- 
quark mass, rrit — 172.5 GeV/c^. 

Any quantity which can be expressed as a function of 
the momenta of the top decay products can be estimated 
as a probability density. We consider the invariant mass 
of the top-antitop pair, and derive the following pdf of 
X, the reconstructed ti invariant mass: 



p{x 



) = f J2^iPk\mt)W{j\pk)S{x -M{pk))dpk 



(2) 



M(pfc) is the invariant mass of the top-antitop pair, given 
parton configuration p^, and W{j\ph) is a transfer func- 
tion which maps the observed jets j to the parton level. 
The sum is over the number of permutations, or jet- 
parton assignments, possible given the observed event. 

The transfer functions W^(j|pfc) used in Eq. ^\ are es- 
timates of the primary parton energy corresponding to 
the observed jet E't- These estimates are derived from 
Monte Carlo simulation, where jets are matched to par- 
tons within Ai? = ^^^ + ^'t> < 0.15. No other jets 
or partons may be within AR < 0.6 of the jet /parton 
match. We separate the response function into 10 GeV 
jet-energy bins and into five bins of jet pseudo-rapidity. 

The matrix-element calculation produces a probability 
density function of the reconstructed ti invariant mass for 
each event. We use this sampled pdf as the observable 
in a likelihood calculation. The probability for an event 
in sample i is a function of the signal fraction fsig in the 
sample: 



P.ifsi 



Jsig-' sig,i + (t Jsigj-' bg,i 



(3) 



Here we regard the hypothetical Z' —>■ ti component as 
signal and SM ti and non-ii processes as background. 
The cross section a for Z' -^ ti and the signal fraction 
are proportional, 






(4) 



where n is the total number of events observed, Ns = 
3 X 2 X 2 = 12 is the number of samples [31], and Li 
and Ai are, respectively, the integrated luminosity and 
absolute signal acceptance. The signal and background 
probabilities Psig.i and Phg,i are, according to their tem- 
plates. 



and 



P,i 



Ph, 



Tsig^^{x)p{x)dx, 



Thg,z{x)p{x)dx. 



(5) 



(6) 
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FIG. 1. Signal templates — summed over &-tag multiplicity, 
jet multiplicity and primary lepton type. 



These templates, T^ig^i and Tbg.i, are the summed event- 
by-event pdfs from each of the model components (e.g. 
the Monte Carlo models of SM ti, W+jets, etc.). The 
templates and the per-event pdfs are normalized to unit 



1 ^= 



p{x)dx = 1 



(7) 



(8) 



where T^^i is the template for component c in sample 
i and Nc is the total number of events used to model 
component c. 

Templates for each model component are constructed 
by summing the per-event probability densities. In Fig.fl] 
we show the templates for Z' — > ti at several values of the 
Z' pole mass, reconstructed according to the procedure 
described above. Note that the mass peak is generally 
prominent. For the samples where the Z' pole mass is 
1000 GeV/c^ and above, fewer Z' particles are produced 
on-shell. The off-shell Z' is produced with mass accord- 
ing to the energy available from the parton density func- 
tions, but decays to a real ti pair. The invariant mass 
spectrum for on and off-shell Z' — >■ ti depends on the 
Z' pole mass; we are able to distinguish the Z' mass, 
given the observed reconstructed ti mass spectrum, even 
in cases where the Z' mass is very large and the Z' is 
largely produced off-shell. 

To model the background components ALPGEN [32] 
version 2.10 with pythia version 6.2.16 [33] parton show- 
ering is used for 14^-|-heavy flavor and M^-|-jets. The shape 
of the ti invariant mass at next to leading order (NLO) is 
softer than the predictions of leading order Monte Carlo 
codes. To correct for this effect, pythia weighted at 
the generator level by mcfm 34J (version 5.8) is used to 
model the SM ti component. The QCD "fake" compo- 
nent (where a jet is mis-reconstructed as an electron) 
is modeled in data, where jets are selected which are 



electron-like in their identification variables. Compared 
to the signal templates (Fig. IT]), we notice excellent sep- 
aration between background and signal shapes. 

We write the likelihood simply as the product of the 
per-event probabilities, 



LicT)=Y[P,ia)-G{iy, 



(9) 



Vj IS nuisance pa- 



where G is a Gaussian distribution, 
rameter j with expectation Vj and uncertainty aj. We 
include among the nuisance parameters each lepton cate- 
gory acceptance and its error, as well as the lepton trigger 
efficiencies. SECVTX tagging acceptance is modeled by 
applying a scale factor to account for differences between 
data and simulation. This scale factor appears as another 
nuisance parameter. We integrate over the nuisance pa- 
rameters to incorporate these systematic effects into the 
likelihood. 

The measurement is affected by uncertainties inherent 
to our model of the background: the amount of initial- 
and final-state radiation (ISR and FSR) and the un- 
certainty on the jet energy scale [30, of the calorime- 
ter contribute most significantly, while systematic uncer- 
tainties due to parton distribution function and color- 
reconnection effects [351 [SS] are negligible. While the 
value of the top quark mass determines the value at which 
the ti invariant mass spectrum rises [10] , the shape of the 
tail of the m^j spectrum is insensitive to variation of the 
top mass within its uncertainty. This analysis is sensitive 
specifically to variations of the tail of the m^j spectrum — 
uncertainty on the value of the top quark mass is not a 
significant source of systematic error. 

We treat systematic uncertainties due to acceptance 
effects as nuisance parameters (cf. Eq. ^), while we 
convolve the likelihood function to include systematic ef- 
fects due to background model uncertainties. 

Background model systematic uncertainties are esti- 
mated by evaluating ensemble tests where the experimen- 
tal circumstances are reproduced and the parameter in 
question (e.g. initial- or final-state radiation, estimates 
of background shape or composition) is varied within its 
error. The resulting distribution of the maximum likeli- 
hood estimate of signal cross section form the basis for 
these estimates of systematic uncertainty. We take the 
difference in the means of the +la and — Icr distribu- 
tions as the error due to the underlying uncertainty at 
each mass value of the hypothetical resonances consid- 
ered. The final tabulation of estimated systematic un- 



certainties for this analysis is given in Table HI 



These errors are assumed to be Gaussian in nature. 
The likelihood of Eq. ^ is convolved with a Gaussian 
distribution whose width is equal to the quadrature sum 
of the individual systematic uncertainties: 



/•oo 

L\a)^ G{cr,a',Aa)L{a')da'. 

Jo 



(10) 



The distribution of reconstructed tt invariant mass ob- 
served in the data is shown in Figure [2] This histogram 



Z' pole 




mass [GeV/c^] 
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450 


0.210 
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950 


0.012 
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0.014 
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0.016 


1200 


0.011 
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TABLE III. Total background model systematic uncertainties 
at each of the mass points considered. These values appear 
as the width of the Gaussian convolved with the likelihood in 



Eq. 1 10 1 
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FIG. 2. The histogram of total probability density for the 
1366 tt candidate events observed in C — 4.8 fb~^. 



shows the integrated probability density in each ti invari- 
ant mass bin. For each event we observe a distribution; 
as a consequence any single event can contribute proba- 
bility density to more than one bin. The data show no 
indication of resonant ti production. Moreover, the dis- 
tribution observed agrees very well with the background 
model over five orders of magnitude. 

At each value of the Z' pole mass considered we cal- 
culate the Bayesian 95% confidence level (C.L.) upper 
limit on the value of the cross section times branching 
ratio for Z' production and decay to ti. This upper limit 
is the value of cross section times branching ratio which 
covers 95% of the area under the convolved likelihood, 
integrating from zero cross section, where a flat prior is 
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3. Expected and observed 95% CL. upper limit for 
a{pp -^ Z') X BR{Z' -^ it) for C = 4.8 fb"^ of integrated 
luminosity as a function of reconstructed ft invariant mass. 
The solid line indicates the observed limit. The dashed line 
indicates the theoretical cross section for a leptophobic Z' |^ . 
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700 
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750 
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800 
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950 
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0.07 


1200 


0.070 


0.08 


1300 


0.104 


0.13 


1400 


0.134 


0.18 


1500 


0.197 


0.24 



TABLE IV. Expected and observed 95% CL. limits on the 
cross section times branching ratio for narrow resonance pro- 
duction. 



assumed. 

Figure [3] shows the expected and observed limits on the 
Z' cross section times branching ratio at each mass point 
considered. These data are also given in tabular form in 
Table [TV] The green band shows the ±lcr expectation for 



the null hypothesis, the yellow band the ±2(t expectation. 
Also included is a curve indicating the theoretical cross 
section for a leptophobic Z' [B]. 

We have searched for narrow resonant states decaying 
to top-antitop pairs in a sample corresponding to an in- 
tegrated luminosity of 4.8 fb"-'^ of CDF II data. Events 
were selected in the lepton -I- > 4 jets topology with 
at least one jet tagged as coming from a h quark. Monte 
Carlo samples were used to model the appearance of nar- 
row resonant states decaying to tt. A matrix-element 
reconstruction method was applied, and for each event 
a probability density function of the reconstructed it in- 
variant mass was sampled. This formed the basis for a 
likelihood fit to extract the cross section times branch- 
ing ratio for narrow resonance production, where non-ti 
background fractions are constrained according to their 
best estimates, and SM and resonant it components may 
vary according to the data. Systematic uncertainties 
such as those arising from the error on the jet energy 
scale, or uncertainties on parton distribution functions 
were incorporated by convolving the likelihood function 
with a Gaussian with width equal to the estimated un- 
certainty on the cross section times branching ratio due 
to the underlying uncertainty. Systematic uncertainties 
arising from acceptance affects such as trigger efhciencies 
or scale factors are treated as nuisance parameters in the 
likelihood function. The benchmark model of a leptopho- 
bic Z' [51 is ruled out at 95% CL. for Z' masses below 
900 GeV/c2. 
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